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We explain the ubiquity and extremely slow evolution of non gaussian out-of-equilibrium distri- 
butions for the Hamiltonian Mean-Field model, by means of traditional kinetic theory. Deriving the 
Fokker- Planck equation for a test particle, one also unambiguously explains and predicts striking 
slow algebraic relaxation of the momenta autocorrelation, previously found in numerical simula- 
tions. Finally, angular anomalous diffusion are predicted for a large class of initial distributions. 
Non Extensive Statistical Mechanics is shown to be unnecessary for the interpretation of these 
phenomena. 
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Modern non-equilibrium statistical mechanics consider 
the study of kinetic behavior of large ensembles of in- 
teracting particles, in an attempt to relate macroscopic 
properties of a system to the microscopic dynamics of its 
constituent. Recently there has been an increasing inter- 
est on long-range interacting systems . The first reason 
is the broad spectrum of applications: self-gravitating 
and Coulomb systems, vortices in two-dimensional fluid 
mechanics 0, wave-particles interaction and trapped 
charged particles Second, unexpected non gaussian 
distributions |4l and non exponential relaxations for au- 
tocorrelations 1, 0, have recently been observed in 
such systems. This was a great surprise, indeed, at vari- 
ance with short-range interacting systems, as all the par- 
ticles contribute to the local field, mean-field predictions 
are usually extremely good and often exact for the cor- 
responding equilibrium statistical mechanics. In accor- 
dance with this simple picture, usual expectations were 
that these systems should exhibit gaussian distribution 
functions and normal exponential relaxations. 

The recent discovery of non gaussian distributions Q 
led to an intense and productive debate on the applicabil- 
ity of usual Boltzmann-Gibbs statistical mechanics to 
long-range interacting systems (sec Tsallis et al in Ref.Q). 
Such non Gaussian distributions have been fitted using 
Tsallis' distributions Q. The striking asymptotic alge- 
braic behaviors for momentum autocorrelations have also 
been fitted using q-cxponcntial functions 0, , derived 
from non extensive statistical mechanics. 

In this paper, revisiting this question, we present for 
the first time analytical results explaining these numeri- 
cal findings, without using concepts of non-extensive sta- 
tistical mechanics. In a recent paper [l^ . we explained 
the ubiquitous occurrence of non equilibrium distribu- 
tions by linking them with stable Vlasov solutions. Here, 
we complete this argument by explaining that the time 
scale for the evolution of such non Gaussian distributions 
has to be larger than the number of particles A^. In an 
other recent paper, we have derived the Fokker Planck 



equation describing a single particle in an equilibrium 
bath 0. Here we generalize this result to out of equi- 
librium bath distributions. More importantly, we also 
report analytical explicit predictions for the momentum 
autocorrelations at large time, explaining its algebraic re- 
laxation. Finally, anomalous diffusion phenomena in the 
asymptotic time limit were found in these systems and 
also highly debated This concept, raised in a great 

variety of different fields such as geophysics, chemical en- 
gineering and disordered media, was analytically under- 
stood in only very few many-degrees-of-freedom systems, 
even if the question of Levy fiights is very popular nowa- 
days. Starting from the microscopic dynamics of a long- 
range interacting system, we present in this letter, ana- 
lytical predictions showing that there exists a full class 
of distributions which leads to anomalous diffusion. We 
give an explicit prediction of the anomalous diffusion ex- 
ponent in terms of the bath distribution function. We in- 
ssit that such an exact prediction of this exponent, from 
the microscopic dynamics of a Hamiltonian N particle 
system is not common at all. 

In this letter, we consider these questions for the 
Hamiltonian Mean-Field model ^3 
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1 
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since it is nowadays thought to be the simplest model to 
study dynamical and thermodynamic properties of sys- 
tem with long-range interactions In addition to its 
pedagogical properties, it corresponds to a simplification 
of one-dimensional gravitational interactions and is an 
excellent first step before the Colson-Bonifacio's model 
for free-electron Lasers • Note that the factor 1 /N is 
the appropriate and classical mean field scaling, relevant 
for long-range interacting systems |l2l H^: indeed, the 
physically interesting limit for such systems amounts to 
let the number of particles go to infinity at fixed volume, 
by contrast with the usual thermodynamic limit. 
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Let us consider the kinetics of the HMF model us- 
ing the Klimontovich approach. We have checked that 
an asymptotical expansion of the BBGKY hierarchy 
leads to the same results. The state of the iV-particles 
system can thus be described by the discrete single 
particle time-dependent density function fd{t,9,p) ^ 
5 {6 — 9j (t)) 5{p~ pj (t)) , where S is the Dirac 
function, the Eulcrian coordinates of the phase 

space and {6i,Pi) the Lagrangian coordinates of the par- 
ticles. The dynamics is thus described by the Klimon- 
tovich's equation 



dfd 
dt 



djd 

de 



dVdfd 

de dp 



= 0, 



(2) 



where the potential V that affects all particles is V{t, 9) = 
— ^da J^^dp cos{d~ a) fd{t, a,p). This description of 
the Hamiltonian dynamics derived from is exact: as 
the distribution is a sum of Dirac functions it contains the 
information on the position and velocity of all particles. 
It is however too precise for usual physical quantities of 
interest but will be a key starting point for the derivation 
of approximate equations, valid in the limit N large and 
describing average quantities. 

When N is large, it is natural to approximate the dis- 
crete density fd by a continuous one f{t,6,p). Con- 
sidering an ensemble of microscopic initial conditions 
close to the same initial macroscopic state, one defines 
the statistical average (fd) = foi^iP): whereas fiuctua- 
tions of probabilistic properties are of order \/\/N . We 
will assume that /o is any stable stationary solution 
of the Vlasov equation. The discrete time-dependent 
density function can thus be rewritten as fd{t,0^p) = 
fo{S,p) + Sf{t,9,p)/\/N, where the fluctuation 6f is of 
zero average. We define similarly the averaged potential 
(V) and its corresponding fluctuations 6V{t, 6) so that 
V{t, 0) = {V) + 5V{t, 9) / VTV. Inserting both expressions 
in Klimontovich's equation (0) and taking the average, 
one obtains 



dfo 
dt 



dfo d{V) dfo 



de 



de dp 



1 /d5Vd5f\ 



N 



The Ihs is the Vlasov equation. The exact kinetic equa- 
tion (0) suggests that for times much smaller than A^, 
and stationary stable solutions /o for the Vlasov equa- 
tion, the rhs term, corresponding to the fluctuation of 
the mean-field potential, can be neglected. This is con- 
firmed, for any finite time, by the Braun and Hepp theo- 
rem MhJ. The innumerable Vlasov stable stationary 
states, far from equilibrium El (see [13 for a first sta- 
bility analysis), then explains the generic occurrence of 
out-of-cquilibrium distributions. These quasistationary 
states do not evolve on time scales much smaller than iV, 
explaining the extremely slow relaxation of the system 
toward the statistical equilibrium. 

Let us now concentrate on stable homogenous distri- 
butions /o(p); which arc stationary since {V) = 0. Sub- 



tracting Eq. from Eq. ^ and using fd = fo+Sf /\/N, 
one gets 



dSf dSf 
dt de 



dSV dfo 
de dp 



dSV dSf / dSV dSf 



de dp 



de dp 



For times much shorter than y/ N, we may drop the rhs 
encompassing quadratic terms in the fluctuations. The 
fluctuating part Sf are then described, by the linearized 
Vlasov equation (this is an other result of the Braun and 
Hepp theorem 0,^3)- This suggest to introduce the 
spatio-temporal Fourier-Laplace transform of Sf and SV. 
This leads to 



SV{uj,k) = - 



7r(4,i + 4 -i) 



" (5/(0, fc,p) 

dp -r-, 7> (4 

i[pk — cu) 



where 



e(w, k) = l + nk {Sk^i + 4 -i) / dp 



+ 00 



dfo 
dp 



[pk — uj) 



(5) 



is the dielectric permittivity. The evolution of the poten- 
tial autocorrelation, can therefore be determined. For ho- 
mogeneous states, by symctry, {5V{uJi,ki)5V{(jJ2,k2)) = 
except if fci = — fc2 = ±1. One gets, after a transitory 
exponential decay, the general result 



{5V{tu±l)5V{t2,Tl)) = \ j^Au: 



-iuj{ti —t2) 



/o(^) 

k(^,i)r 



(6) 

This is an exact result, no approximation has yet been 
done. 

A similar, although longer, calculation shows that the 
rhs of Eq. identically vanishes at order This 
proves that Vlasov stable distribution function will not 
evolve on time scales smaller or equal to iV, in agree- 
ment with the N^-"^ scaling law which was numerically 
reported This is the first important result: 

generic out of equilibrium distributions evolve on time 
scales much larger than N. 

Let us now consider relaxation properties of a test- 
particle, indexed by 1, surrounded by a background sys- 
tem of (A^— 1) particles with a homogeneous distribution. 
The fluctuation of the potential is thus 

/•2tt /• + 00 1 

SV{t,e) = - da dp cos{e-a)Sfit,a,p) ^ cos (61 - 6ii) . 

(7) 

Using the equations of motion of the test particle and 
omitting the index 1 for the sake of simplicity, one ob- 
tains p{t) = p(0) - J^du{dSV{u,e{u)))/{de)/VN. By 
introducing iteratively the expression of 9 in the rhs and 
expanding the derivative of the potential, one gets the 
result at order 1/7V. The key point is that this approach 
does not use the usual ballistic approximation. As a con- 
sequence, we obtain an exact result at order 1/N. This 
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is of paramount importance here to treat accurately the 
collective effects. As the changes in the impulsion are 
small (of order 1/^/N), the description of the impulsion 
stochastic process by a Fokker-Planck equation is valid. 
This last equation is then characterized by the time be- 
havior of the first two moments {{p{t) — p(0))"). Using 
the generalization of formula © when the effect of the 
test particle is taken into account, one obtains in the 
large t-limit 



{{p{t)-pm) 

{ipit)-piO)f) 



dD 
dp 



^+00 t 
N 



'^-'4 Dip), 



(P) 



]_dfo 
fo dp 



Dip) 



(8) 
(9) 



where the diffusion coefficient Dip) can be written as 



D(p) = 2Re / dte'"' {<51/(t, 1)<5V(0, -1)) 



(10) 

These results are the exact leading order terms in an ex- 
pansion where 1/N is the small parameter. We obtained 
previously equivalent results, but restricted to equilib- 
rium (gaussian) /o 1 by a rather different and instruc- 
tive approach. In this last paper, a comparison of the 
diffusion coefficient Dip) with N particle numerical sim- 
ulations is presented, illustrating that both results are 
undistinguishable. Let us carry on by explicitly evalu- 
ating the diffusion coefficient for a homogenous Gaus- 
sian distribution function /g(0,p) = (3/i2Tr)^e~^P 1"^ . 
In that case, after straightforward calculations, one gets 
the expression derived in 0. The diffusion coefficient 
Dip) has the asymptotic expression y^TT[3/2 e~^P 1"^ for 
large p values (see Fig. 

Interestingly, the method presented here, can be used 
for any Vlasov-stable out-of- equilibrium distributions. 
For instance, in Fig. ^ wc present the result for the wa- 
terbag distribution. 

Using time variable r = t/N as suggested by Eqs. © 
and the Fokker-Planck equation describing the time 
evolution of the distribution of the test particle is 



N f 9fiiT,p) 1 <9/o , . . 



dfiiT,p) _ d 
dr dp 

(11) 

We stress that this equation depends on the bath dis- 
tribution /o. It is valid both for equilibrium and out-of- 
equilibrium /o , provided that fo is a stable stationary so- 
lution of the Vlasov equation. In the limit r ^ 00 (more 
precisely 1 ^ r <C N), the bracket vanishes: the pdf /i 
of the test particle converges toward the quasi-stationary 
distribution /o of the surrounding bath. This is in com- 
plete agreement with the result that /o is stationary for 
times scales of order N. 

For a large class of quasi-stationary distributions fo, 
using Eq. Hll|) . we can compute the momenta autocorre- 
lation (p(t)p(O)), fitted numerically HQ, [ill with power 




FIG. 1: Diffusion coefficient -D(p) in the case Hn/N — 2 
for a Boltzmann tliermal bath (solid line) and a waterbag 
distribution (dashed line). 
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FIG. 2: The solid line represents the time evolution of the 
momentum autocorrelations obtained using the numerical in- 
tegration of the Fokker-Planck equation for a test particle in 
a surrounding gaussian reservoir (/3 = 0.25). Its slope can 
hardly be distinguished from a 1/r law (dotted hne). 



laws, stretchcd-exponcntials or q-cxponential. The sec- 
ond important result of this letter is that the time de- 
pendence of the momentum autocorrelation function scale 
with N as (p(t)p(O)) = Cit/N) where C(.) is a function. 
Let us first present the particular but very important 
case of a test particle in contact with a gaussian distri- 
bution fo (equilibrium bath) . Fig. |2] shows the momenta 
autocorrelation (|3(t)p(0)), numerically computed from 
the Fokker-Planck Eq. (|11|) : it presents an unexpected 
very slow relaxation which can, numerically, hardly be 
distinguished from a l/r-law. As shown below, the very 
fast decrease of the diffusion coefficient shown in Fig. 
is actually the key point in these interesting and unusual 
properties of the momenta autocorrelation. 

One derive analytically the large time behavior of the 
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autocorrelation function, for any Vlasov stable distribu- 
tion /o, u sing an asymptotic expansion inspired by 
(see also |2y|). By introducing the appropriate change 
of variable x = x(p), defined by Ax /dp = D(jp), one 
ends up with the constant diffusion coefficient Fokker- 
Planck equation 



dh{T, x) 
dr 



d d.h{T,x) , a^y, \ 



In (/Db)/o(p)) . Using e{p, 1) 1, 



where ^'(2;) 

Eq. Iini) shows that D{p) n'^foip). 

Let us first consider distribution functions so that 
foip) '^'^ C* oxp(— 7^"^), which includes not only the 
gaussian (S = 2) and exponential tails (d = 1), but also 
stretched-exponential ones with S > 0. Asymptotic anal- 



ysis leads to p{x) 



ip{x) a In \x 



2/7 (In(x))^^'' whereas 

with a ~ S. (13) 



Interestingly, this log-potential is a limiting case between 
a pure diffusive process leading to self-similar solutions 
and exponentially decreasing /i-solutions in a strongly 
confining potential. Asymptotic behavior is still 
valid for distribution functions with algebraic tails, i.e. 

fo{p) ^^^r^°° Cp~^ where v > i. However, a ~ iv/{2 + v) 

&ndp{x) 

By considering the ansatz fi{x,t) = (px{x)e~^'^ , one 
gets that the spectrum of the Fokker-Planck operator 
corresponds to only one bound state fi{x) ~ /o, normal- 
izablc for a > 1, with an associated eigenvalue located at 
the bottom of the continuum: the absence of gap forbids 
a priori any exponential relaxation. 

As we arc interested in the asymptotic largc-r limit, 
we will restrict the analysis to the small-A regime. A 
matched asymptotic procedure in the two distinct re- 
gions, |a;| > £ and |a;| < £, solves the equation. In the 
first domain, introducing z = ^/Xx and gx{z) = zip\, 
and using the asymptotic result (|13|l . one ends up with 
z'^g'x + a'x{a - 2)z + gx{2 - 2a + z^) = 0. The solu- 
tions can be expressed in terms of Bessel functions of or- 
der ly, and Y^, as gx{z) ^ ^A,f 2:^^""'^^ J(q+i)/2(^) + 
Bx^e z^'^~°'^^^Y{a+i)/2iz)- In the domain \x\ < £, where 
one neglects the term proportional to the vanishing 
eigenvalue A, the solution is (px{x) ~ Dx.ee 
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x,i e 



-ip{x) 



In order to compute the mo- 



menta autocorrelation function, noting that p(x) is an 
odd function of x, we focus on odd eigenstates, obtained 
by considering Dxj = 0. By taking care of the match- 
ing condition in .t = £ and of the normalization con- 
dition, one ends up with the scaling: Ax^i %/A/2, 

Bx,t A/4 and Cx,t C(0A(5/2). 

All these results are finally useful to derive the mo- 
mentum auto-correlation. Indeed using as initial con- 



dition /i(x,0) = N{£)p{x)e-^^-') = d\p(X)^x[x) 
where N{X) is the ground state normalisation factor and 
^(A) = N{£) J dxp{x)(px{x), one gets 



(p(r)p(O)) 



1 



dA^(A)'e 



2„-At 



(14) 



Ni£) 

Its limiting behavior in the r — > 00 limit will be given by 
the behavior of /x(A) when A ^ 0, which is itself deter- 
mined by the large behavior ofp{x). 

Let us be more specific in several important cases. The 
choice p{x) oc (Inx)^/'' which corresponds to distribu- 
tion functions with gaussian, exponential or stretched- 
exponential tails, leads to the important result 



(p(r)p(O)) 



-.+cc^ (lnT)2/« 
cx . 



(15) 



independent on a. This very slow algebraic relaxation, 
with logarithmic corrections, agrees with numerical sim- 
ulation shown in Fig. [3 for the Gaussian case 5 = 2. 

An extension of this approach to distribution function 
fa{p) with algebraic tails is of prime interest. Following 
the same steps in the case p{x) ^ x^ {rj < {a — l)/2) 
leads to /i(A) ^oc^ ;\("-2'?-3)/4 ^^^^ ^^-^^ limiting be- 
havior (a;''(T)x''(0)) ^lx°° r''-("-i)/2. Applying this re- 
sult to the algebraic decay of the distribution function 

fo{p) ''''^ Cp^^ , with u > 2 and using the correspond- 
ing asymptotic expression for p{x), discussed above, we 
get = 2/{2 + v) and a = ivji^ ^v). We thus obtain 
the algebraic decay 



(p(t)p(O)) '"Ix °° t' 



(16) 



This is the third important result: the equilibrium distri- 
bution and generic out of equilibrium distributions lead 
to algebraic large time behaviors of the momentum auto- 
correlation functions. 

Finally, from the momenta correlations, one usually de- 
rives the angle diffusion {{0{t) — d{Q))'^) = 2De r where 
Dg is defined via Kubo formula Dg = J^°° dr (p(t)p(O)). 
However, asymptotic result (|15|l shows that this inte- 
gral diverges leading to small anomalous diffusion for 
stretched exponential bath distributions. By contrast, 
for distributions with algebraic tails, since the exponent 
(i' — 3)/(2 -I- 1/) in Eq. (|15|l is smaller than one, we predict 
strong anomalous diffusion. This is the fourth important 
result. 

In summary, this work confirms and explains the un- 
expected kinetic behavior of long-range interacting sys- 
tems: ubiquity of non gaussian distribution and algebraic 
behavior for momenta autocorrelation functions. Using 
traditional kinetic theory, we also predict strong anoma- 
lous diffusion for angles for a large class of initial distri- 
butions. Finally, let us anticipate that several physical 
systems with long-range interactions should exhibit sim- 
ilar features. One would in particular quote dynamics of 
vortices and cold atoms physics 
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